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ABSTRACT 


The present investigation deals with flexural-torsional 
buckling of simply supported, two way linearly tapering, 
monosymmet ri c I beam-columns/tie— beams under four types of loading 
(i) Moment gradient, (ii> Uniformly distributed load, (iii> Single 
point load, (iv) Two symmetrically placed point loads. In the 
case of transverse loading, the loads are assumed to be applied at 
the top flange, at the shear centre or at the bottom flange. 
Results are presented in either graphical or tabular form. It is 
shown that the design specifications of 13:800-1984 could be 
unsafe in predicting the buckling capacity of monosymmet r i c 1 


beams in some cases. 
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CHAPTER I 


INTRODUCTION 

1.1 STABILITY 

The stability of a structure essentially means the stability 
of its equilibrium configuration or state. In a practical sense, 
an equilibrium state of a structure or a system is said to be 
stable if accidental forces, shocks, vibrations, eccentricities, 
imperfections, inhomogeni ties or other probable irregularities do 
not cause the system to depart excessively or disastrously from 
that state. In a mathematical sense, the stability is usually 
interpreted to mean that infinitesimal disturbances will cause 
only infinitesimal departure from the given equilibrium 
configuration . 

The loss of stability under compressive loads is usually 
termed as structural <or geometric or form) instability, commonly 
known as buckling. There are three basic ways by which 
instability can occur in a single structural member: 

- Pure flexural buckling about the major principal axis. 

- Pure flexural buckling about the minor principal axis- 

- Pure torsional buckling about the shear centre axis. 

Depending on the sectional properties and external loads, a. 

structural member becomes unstable in any of the above three 
individual buckling modes or by a combination of these buckling 
modes. For example, consider a structural member with 



monosymmetric I cross section subjected to axial and transverse 
loads only. The instability can occur by 

— Pure flexural buckling about the major principal axis 
(in-plane buckling). 

— Combination of flexural buckling about the minor 
principal axis and torsional buckling about shear centre 
axis (flexural-torsional buckling). 

The buckling capacity of such a member is the elastic 
critical load obtained from analysis of above two modes. 

1.2 FLEXURAL T<»£TONAL BUCUXING 

When an I beam is subjected to axial and transverse loads, it 
may buckle sideways if its compression flange is not laterally 
supported. It happens because the compression flange, which is in 
effect like a column on an elastic foundation, becomes unstable. 
At the critical loading there is a tendency for the compression 
flange to bend sideways and for the remainder of the cross section 
to restrain it from doing so. The net effect is that the entire 
section rotates and moves laterally. Thus, lateral buckling is a 
combination of twisting and lateral bending brought about by the 
instability of the compression flange. Therefore, it is termed as 
f lexural— torsional buckling. Since thin walled open sections are 
most susceptible to such type of buckling, the discussion in 
subsequent articles will be limited to these shapes, with 
particular reference to monosymmetr ic I sections. 

The analysis for flexural-torsional buckling is usually based 
on following assumptions: 

- The material is homogeneous, isotropic and elastic. 
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- Plane sections originally perpendicular to he 

longitudinal axis of the member, remain plane and 
perpendicular after bending. 

— Flexural deformations are small. 

The above three assumptions enable the use of classical beam 
theory for the purpose of analysis. 

— The cross section of the member is allowed to rotate and 

translate in a plane perpendicular to the longitudinal 
axis of the member, but no distortions in cross section 
shape is allowed. (This assumption ignores the effect 
of local buckling of flanges and web.) 

— The beam is initially undistorted and unstressed, (i.e. 

no geometric imperfactions and no residual stresses.) 

Within the limit of above assumptions the buckling capacity 
of a member depends on the following factors! 

— Section properties. 

— Loading conditions. 

— Lateral restraint provided and its position along the 

span. 

— Displacement, rotational, torsional and warping 

constraint at the ends. 

— Variation in shape of cross-section along the span. 

1.3 LITERATURE REVIEW 

The elastic flexural-torsional buckling of uniform doubly 
symmetric members has been studied at length. Numerical solutions 
and approximate buckling formulae for such members, under various 
loading and restraint conditions, are available in standard texts 



e.g. Timoshenko and Gere (1961), Galambos (1968), Chajes (1974), 
Trahair and Bradford (1988). 

Much of the research work on the buckling of tapered members 
is focused on doubly symmetric web tapering members. SoiTie early 
work includes inplane buckling of 1 sections under axial thrust 
and axial thrust-end moments combinations with various boundary 
conditions (Fogel et al 1962; Gore et al 1962 Gi ri javallabhan 
1969). The behaviour of linear and parabolic web tapering 
cantilever beams under external torque was studied by Lee and 
Szabo (1967). Culver and Preg (1968) derived differential 
equation for flexural-torsional buckling of one way web tapering 
beams and solved it using finite difference method for moment 
gradient load case. Kitipornchai and trahair (1972, 1975) studied 
the buckling behaviour of two way tapering members under a central 
point load. This study also includes effect of variation in the 
thickness/width. These results were confirmed by Bradford and Cuk 
(1988) who investigated lateral buckling of monosymmetr i c I beams 
using the finite element method. Brown (1981) studied 
lateral-torsional buckling of two way web tapered I beams. He 
derived the governing differential equation using energy approach 
and solved it using finite difference method for central point 
load. 

One of the difficulties encountered in the design of 
monosymmetr i c I beams has been the determination of section 
properties required for calculating the elastic critical moment. 
Kitipornchai and Trahair (1980) overcame this difficulty by 
introducing some new parameters. Using these parameters and the 
energy approach, Kitipornchai et al (19S6a, 1986b, 19S8a) studied 



the buckling of uniform members for various load cases. To 
simplify the flexural-torsional analysis of beam-columns / 
tie-beams, Wang and Kitipornchai (1989) proposed a few additional 
parameters and studied the behaviour of uniform mono symmetric 
beam-columns/tie-beams under various loading conditions 
(Kitipornchai and Wang, 1988b). 


1.4 AIM OF IMVESTIGATION 

15:800 (1984) uses the Merchant Rankine formula to get the 
maximum permissible bending compressive stress designing 
a member: 


& 


0.66 f . f 
cb y 


be 


C(f . )" + (f 
cb y 


where f^ is the yield stress of the steel; n is the imperfection 
factor and f . is the elastic critical stress in bending. To 
determine f^^ the code proposes a formula or suggests the use of 
an elastic flexural-torsional buckling analysis. The formula 
proposed in the code is derived by making use of the equation 
which is obtained for the elastic critical moment for a simply 
supported, uniform, monosymmetric, 1 beam subjected to a constant 
moment. In actual practice the member can be nonprismatic and the 
loading can be of any type resulting in various forms of bending 
moment variations. Under this circumstances one is not sure 
whether it would be correct to use the formula given in 15:800. 

Two way linearly tapering members, which give effective 
utilization of structural materials, are widely used for plate 
girders and gantry girders. In such members, only the depth 


varies along the span while the area of each flange remains 
constant. With this in view the aim of present investigation is: 



- to study the flexural-t 
members under an axial fo 
external loads: 

* Linear moment gradien 

* Uniformly distributed 

* Single point load pla 

» Two symmetrically pla 

- to give values of nondimen 

by means of which f . , and 

cb 

sections available in prac 

l.B OUTLINE OF THESIS 

The thesis consists of three chapters. Chapter 1 is the 
introduction. In chapter S, all the section properties have been 
defined. The formulation of the buckling problem for various 
loading conditions is also given in this chapter. Results of the 
analysis are presented and discussed in Chapter 3. Whenever 
possible results obtained from the present analysis are compared 
with the provisions in 15:800. Tables are presented for doubly 
symmetric I sections and various loading conditions which can be 

cb* 


orsional behaviour of such 
rce combined with following 

t. 

load. 

ced anywhere on the span, 
ced point loads. 

sional elastic critical moment 

thus ct , can be obtained for 
be 

t i ce . 


used for determination of elastic critical stress f 



CHAPTER 11 


PROBLEM FORMULATION 

2.1 I NTRODUCTI ON 

In this chapter, a general formulation is presented for the 
flexural— torsional buckling analysis of two way linearly tapering 
monosymmetr i c I beam-columns/t ie-beams . The method of 

investigation is similar to that adopted in the recent work of 
Kritiponchai and Wang (1988b). The various load cases which have 
been considered in this investigation are shown in Fig. 2.1. 

To simplify the analysis, following assumptions are made in 
addition to those listed in Sec. 1.2: 

- Resultant of transverse loads act solely in the vertical 
plane passing through the shear centre of the section. 

— Resultant of axial loads lies in the plane of web. It 
is parallel to the top flange and acts through the 
centroid of the section at supports. 

— No intermediate lateral restraint is provided. 

- The displacements and moments about major and minor axis 
of the section vanish at the ends of the member. The 
torsional boundary conditions are, zero twist and zero 
warping restraint at the ends of the member. 

- The flange width, the flange thickness and the web 
thickness are constant along the span, while the depth 
varies linearly with maximum value at the centre. 

- The angle of tapering is small. 
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A typical structural member with general monosymmet ic 
1-section is shown in Fig. 2.E. The length of the member is EL. 
The distance between the flange centroids is ah at the «nd5 and h 
at the centre. The width and thickness of the top and bottom 

flanges are designated by T^ and Tj^, respectively. The 

origin of the coordinate system is at the left support coinciding 
with the centroid of the section with z axis along the length of 
the member and the positive x and the y axis as shown in the 
figure. The member is assumed simply supported at the ends and 
Mv\,restrained against warping, i.e. it can freely rotate about the x 
and y axis at the ends but prevented from rotating about the z 
axis at these points. 

The Rayleigh-Ritz energy method has been used in the 

analysis. It involves minimizing the total potential energy 
(which is the sum of the strain energy U, stored in the buckled 
member and the potential energy V of the external loads) with 
respect to the buckled shape. We can write the following 
expression for the strain energy stored in the member when it 

buckles out of its plane. 


U 


= E 

i=1,E 




El 


yzi 


dz^ 


dz 


is 


GJ 


zi 


dz 


dz 


+ 4 r El . dz "I 

E J wzi . Z J 

dz 


(E. 1) 


In the above expression u is the lateral displacemenht and <jlr is 
the angle of twist. The first term on the right hand side 
represents the strain energy due to bending in the lateral 
direction, the second term is the contribution due to torsion and 
the last one is due to warping. The suffix i indicates the 


Dortion of member for which the expression is valid. 


Note that 



L 



L 


Fig. 2.2 (a): Typical web tapering structural member. 



Fig. 2.2(b) Monosymmetric I -section. 



for i = 1, integration limit is 0-L, while for i = £, integration 
limit is L-EL. 

The expression for the potential energy V for different load 
cases is given in Sec. 2.3 and 2.4. 

2.2 SECTION PROPERTIES 

The determination of the section properties required for 
calculating the elastic critical moment of a monosymmet r i c section 
is not straight forward, and the effort required is prohibitive in 
routine design. Kitipornchai and Trahair (1980) gave following 
four parameters for prismatic members. These parameters help in 
simplifying the analysis. 

I . I . 

(i) p = -= ^ O < p < 1.0 (2.2) 

■•■yt \b ■*^yo 


where I . and I . are the second moments of area about the y axis, 
yt yb 

of the top and the bottom flanges, respectively. 

Apt A. 


(ii) p 


A 4- A 
"Ft "Fb 


Ft 


0 < p < 1.0 


(2.3) 


where Ap^ and Apj^^ are the areas of the top and bottom flanges, 
respectively. The parameters p and p can =be viewed as a measure 
of the degree of section monosymmetry. Note that p = p = O 
represents an inverted T— section, while p = p = 1.0 denotes a 
T-section. For a doubly symmetric I section, p = p = 0.5. 

(iii) When a monosymmetric section subjected to transverse loads 
twists during buckling, the major axis moment M causes an 
additional torque which results in an effective change in the 
member torsional stiffness from GJ to (GJ + M in which is 
the monosymmetry property; the expression for is! 



!j! = 0.9 <EP- 1 > 

h V X 

or. = 0.9 < 2 p - •>’ '‘y ^ ^ 

fL t 90 way taperins n.a»>.ars. tho atova exprasaion can ba 

slightly modified as! 


ft . = 0.9 ( 2 p- 1 > ^*zi " - 
' xzi 


( 2 . 4 ) 
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16 GJ^L 


( 2 . 6 ) 


, is defined as the ratio of web area to flange 

(v) Parameter 

+ A- fiance z from the origin. It varies linearly along 

area at any distance 

j .w. -L 1 % i3> w r i 1 1 e n as 
the span and can oe wi 
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z1 
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where, 
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1 + ( 
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T. 01 
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Using above parameters, the following useful expressions 
derived: 

(vi) Minor axis moment of inertia 

f r) . T ^ ^ 

I . = I f '• + ^ p I 

yo I + (1-^) ' 

where , 

I = I . + I . . 
yo yt yb 

For thin walled members, the second term in the parenthesis 
can be safely neglected. Then, 


>}] 

can be 


I . = I 
yz 1 yo 

(vii>Torsion constant 


J . 
zi 


1 + 


"zi ■'w 
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where , 


J = 
o 


®t . »b 


3 3 

(viii) Polar radius of gyration about the centroid 
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<ix) Distance between centroidal axis and shear centre axis 


f ^ _ I ^ 

^ozi [ 2 (1 j ^ J "zi 

(x) Eccentricity of centroidal axis 

= { 

(xi) Warping section constant 
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wzi yo zi 


( 2 . 11 ) 


r 2 " "'zi 1 

t h - i 

- 2 ( 1-/J ) + rri ^ 

1 2 J 

r ^zi 1 

[ 2 (Hrr)) J 


( 2 . 12 ) 


(2.13) 


2.3 MOMENT GRADIENT LOAD CASE 


Consider a two way web tapering structural member subjected 
to moment M at one end and /9M at the other end as shown in Fig. 
2.1(a) . 

The expression for the potential energy V of the inplane 
moment and the axial force when the member buckles out of its plane 
is given by. 


V = r -f r M . ( ) ip dz +ir M .ft . ( -^ )^ dz 

. — „ 1 J xzi . 2 ^ zJ xzi ‘ xzi dz 

1=1,2 '• dz 

[ 


, ^ . ,2. < ^ ,2 - Ey . 

dz zi dz "^ozi 


( ^ > ^ ] dz } 


( 2 . 14 ) 

As indicated earlier, the suffix i denotes the portion of inember 
for which the expression is valid; for i = 1 integration limit is 
0— L, while for i = 2 it is L— EL. In the expression (2.14), we 
have , 

2 

= Polar radius of gyration about shear centre 



(E.15) 
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= Moment about centroid of the section at any 
distance z from the origin 
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It is found convenient to use the following nondimensional 
terms s 
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The total potential energy FI = U + V may be non 
dimensionalized by denoting 

n = EOL / El , therefore, 

yo 

n = (U + V) (2.18) 
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Using Eq . (2.1) to (2.17) in Eq. (2.18), the nondimensional 
expression for total potential energy can be written ass 
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As mentioned in Sec. E.1, the boundary conditions ares 


d u 

y — O at z — Ot E.O 

dz^ 

i.e. the displacement and the moment vanish at the ends of the 
member . 

^ = O at z = O, 2.0 

dz 

which imply zero twist and zero warping restraint at the ends of 
the member. 

In view of the boundary conditions mentioned above, the 
nondimensional displacement and rotation functions, u and , may 
be parameterized by the following trigonometric series. 

_ n — _ n — 

u = E sin = E e sin (2.20) 
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The Rayleigh— Ritz method requires 
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Substituting the relation (2.20) in the expression (2.19) and 
performing above operations, we get the following two Ritz 
equations . 
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n matrix form, the preceding two equations can be written ass 
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where, C A □ is the stiffness matrix and C B D is the load matrix. 
The elements of matrices CAD and C B D are! 
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where, i = 1, ..., n; j = 1, n. (E.28) 

A nontrival solution of Eq. (2.E3) exists when, 

|CA3-kCBD| = 0 

This yields a power series of order (n + n) whose lowest possible 
root is the critical value of the buckling parameter k. 
Sufficient number of terms (n and n) in the trigonometric series 
(Eq. E.EO) must be taken for accurate results. 



2.4 TRANSVERSE LOAD CASES 


Following the procedure outlined in the previous section, 
elements of the stiffness matrix C A □ and the load matrix C B 3 
can be obtained for load conditions shown in Figs. E.1 (b,c,d). 

In each load cases, following load positions are ana;\sed: 

- Top flange loading 

- Shear centre loading 

— Bottom flange loading 

At this stage it is necessary to define a parameter ^^i' 

which describes the load height below shear centre. The values 

of a . are listed below: 
zi 


a 


zi 


-(1-p) h^^ ; top flange loading (TF) 

0 ; shear centre loading (SC) (E.29) 

joh^^ ; bottom flange loading (BF) 


The elements of the matrix CAD depend on the geometry of the 
structural member, material properties and the axial load. 
Therefore, the matrix CAD will remain the same for all 
transverse load cases. On the other hand, the matrix C B D will 
change with loading conditions. 

Instead of writing each step, only those expressions, which 
change with loading conditions/load positions are given below. 
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It may be noted that a . /h .) is a constant and its value is given 

Z 1 SI 

by Eq. (2.28). 


2.4.2 Single Point. Load 
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2.4.3 Two Point Load 
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A computer program is developed to -formulate the matrix C A 3 

and the matrix C B 3. Numerical integration scheme is adopted for 

evaluating h. .and f. . elements of the matrices. NAG Subroutine 
iJ xj 

FOEbjf is used to solve the eigenvalue problem of the form 
CA3{;X3 = kCB3tX>. 

The results of the analysis are presented and discussed in 


the next chapter 



CHAPTER III 


RESULTS AND DISCUSSIONS 


3.1 I NTRODUCTI ON 

In Chapter II, expressions for elements of the load matrix 
and the stiffness matrix were derived for various load cases. As 
seen from Eq. (2.27), the expression for f^ elements involves , 
Bj^, h and some combinations of these parameters. It is 
more convenient to use the clear height of web h , instead of h 
which is the distance between the flange centroids. Using h^, the 
web/flange area ratio X) at the centre can be written as 


(T /h ) 
w w 

<B./h ) <T*/h ) 
t w t w 


(1-^) (T /h ) 
w w 

(B./h ) <T./h ) 
b w b w 


(3.1) 


If u, p, (T /h >, (B/h ). (T/h > are known, then the section 
w w w w 

dimensions can be uniquely determined for a given h^. Here B and 
T are the width and the thickness of the governing flange, i.e. 
the flange having higher moment of inertia about the minor axis of 
the member. (Top flange for p > 0.5 and bottom flange for p<0.5). 
For p = 0.5, flange having greater area will be the governing 
flange . 


In 

general , 

the stability condition for 

beam- 

“Columns 

/ 

tie-beams should 

be influenced by the 

parameters p. 


K y CJI 

and 

(T /h ), 
w w 

(B/h ), 
w 

(T/h ). However, 

w 

in practice 

the 

range 

of 


variation in the values of last three parameters is very small. 

For example, for ISMB sections the value of u is 0.5? (T /h ) is 

w w 



in the range 0.022 to 0.032; (B/h^) varies from 0.38 to 0.56 and 

<T/h ) varies from 0.04 to 0.06. Therefore, it will be justified 
w 

if we keep p , (T /h ) , CB/h ), (T/h ) constant and study the 

w w w w 

flexural— torsional buckling for various combinations of p, a, K 
and k only. For the purpose of present analysis, following values 
are adopted: 

pi = 0.5, (T /h ) = 0.025, (B/h ) = 0.467, (T/h ) = 0.04. 
w w w w 

Beam-columns/tie-beams with cx = 0.25, 0.50, 0.75, 1.0 and p =0.1, 

0.3, 0.5, 0.7, 0.9 have been studied for various load combinations 

with axial load ratio varying from -1.0 to +1.0. The value of the 

parameter K varies from 0.0 to 2.0. 

Convergence studies were conducted to determine the number of 

terms in the trigonometric series (Eq. 2.20) required for accurate 

solutions. It was found that for prismatic members with 

syymmetric loads, n = n = 3 terms are sufficient. However, for 

tapering members and unsymmetric loading, atleast 7 terms are 

required. The numerical results presented herein are based on 

n = n = 7. 

3.2 MOMENT GRADIENT LOADING CASE 

A typical stability criteria curve for monosymmetr i c, 
prismatic beam-columns with p > 0.5 and subjected to uniform 

sagging moment is shown in Fig. 3.1. The point A corresponds to 
flexural— torsional buckling load for columns. If the combination 
of X and k is such that it falls inside the area enclosed by the 
curve, the member is stable. When the axial load ratio exceeds 
the value X^, member is unstable unless accompanied by adequate 
applied moment within the stability criterion (i.e. between points 
B and C in Fig. 3.1). The applied moment serves to stabilize the 




Fig. 3.1 Typical stability criteria curve for 
monosymmetric beam-columns. 


in the 


member by decreasing the compressive normal stresses 
smaller flange which is further away from the shear centre; 
thereby decreasing the destabilizing torque caused by the 
compressive stress components as the section warps during buckling 
(This is the Wagner effect). For stability, the applied mioment 
must be sufficiently large such that the resulting stresses in the 
smaller flange become tensile so as to provide an effective 
restoring torque (Fig. 3.2). On the other hand, when the applied 
moment is too large, the member will become unstable because the 
larger flange is not able to sustain the increased compressive 
stresses (above point C in Fig. 3.1). The same explanation is 
valid for prismatic beam-columns with p < 0.5 and subjected to 

uniform hogging moment. Results for moment gradient loading are 
presented in Figs. 3.3 to 3.10. 

Figures 3.3 to 3.6 give stability criteria for a typical 
monosymmetric tapering member with K = 1.0, ot = 0.25, 1.0 and p = 
0.1, 0.5, 0.9: values of ft considered are -1.0, 0.0 and 0.5; axial 
load X varies from -1.0 to +1.0. Figures 3.7 - 3.9 show influence 
of a and p on the buckling capacity for K = 1.0, X = -1.0, 0.0, 

0.5 and ft = -1.0, 0.0, +1.0. The Fig. 3.10 shows influence of K 
and p on buckling capacity for oi = 0.25, 1.0, X = 0.25, 0.75 and 

for ft = 0.5. From these results following points can be noted: 

(i) For single curvature bending of prismatic members (Figs. 3.3a 
and 3.4a), the buckling capacity k increases with increasing 
p value for the sagging moment case, while for the hogging 
moment k increases with decreasing p value. This may be 
attributed to the Wagner effect which is more favourable when 
the flange with the smaller is in tension during buc ling 
(Kit i pornchai and Wang, 1988b). 
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Combined axial force 



provide effective restoring torque 

Fig. 3.2 Influence of tensile stresses in smaller flange. 



(ii) The trend similar to that in <i) is observed for tapered 

members when the axial load is compressive (Figs. 3.3b and 

3.4b); however, it is not so when the axial load is tensile 

and sufficiently large. the axial tension P causes the major 

axis moment This component acts in the same direction 

as the externally applied moment and therefore has a tendency 

to destabilize the member. When the axial tension is small 

this component is over shadowed by the Wagner effect, but it 

becomes significant when the axial tension is large. That is 

why the buckling curve for p = 0.5 is higher than that for p 

= 0.9, when axial tension is very large. The effect is even 

more severe for p = 0.1 and for hogging moment case, where in 

addition to Py , Pe also acts in the direction of the 
02 z 

externally applied moment. 

(iiilFor double curvature bending ifi > 0), the Wagner effect has a 
mixed influence, since smaller flange at one end is in 
compression while at the other end it is in tension. Figs. 
3.5a and 3.6a show that buckling capacity for prismatic 
member will be greatly curtailed if axial compression is very 
high. The buckling behaviour is governed by the portion of 
smaller flange which is under compression. As expected, for 
ft = +1.0 the buckling capacity for p and <1-p) is the same. 
The buckling capacity is greatest for doubly symmetric 
sections (Fig. 3.6a). 

(iv) Figures 3.5b and 3.6b show stability criteria for tapered 
members under moment gradients with ft — +0.5 and ft — +1.0, 
respectively. For p = 0.1 and high axial compression, a 
significant reduction in the buckling capacity takes place. 

But it is not so for p = 0.9 because the component Pe of the 

2 *. 
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major axis moment reduces the portion of smaller flange which 
is subjected to compressive stresses. For ft = +1.0 and with 
axial tension (-1.0 < X < 0), the buckling capacity for p 
and (1-p) is almost the same. 

(v) Figure 3.7 shows behaviour of tapering beam (X = 0) under 

moment gradient case. For single curvature bending ift < O), 

as the value of a increases, the buckling capacity reduces 
/otO valiL&s f ajY\d /'ncveassw 

for^high values of p. This is again due to Wagner effect. 
For p < 0.5, the smaller flange is under compression which 
has a tendency to destabilize the member by reducing 
effective torsional rigidity. As the value of a increases, 
the distance of smaller flange from the shear centre 
increases which further increases the destabilizing torque. 
Thus buckling capacity reduces with increasing a. Just the 
reverse is the case when p > 0.5. 

(vi) For double curvature bending of beams {ft — +1.0, Fig. 3.7), 
the buckling capacity for p and (1— p) is equal and is highest 
for a symmetric section. It is interesting to note that for 
extreme values of p, the Wagner effect is more prominent. As 


the section monosymmetry reduces (p approaches 0.5), Wagner 
effect diminishes. 

(vii)For single curvature bending of beam-columns {ft < O, Fig. 

3.8), the buckling capacity decreases with decreasing a. for 

all values of p. for p < 0.5, Py and Pe components of the 

0 2 2 

major axis moment acts in the same direction as the 
externally applied moment. As ex decreases, Pe^ component 
increases, which increases the effective moment at the centre 



of the member. Thus buckling capacity decreases. For p > 0.5, 

Py component of major axis moment reduces the eternally 

O Z ■ 

applied moment. As a decreases, the effective moment near 
the end of the member (where cross section is least) 
increases. Thus the buckling capacity reduces. 

(viii) The behaviour of beam-columns under double curvature 
bending {ft = +1.0 Fig. 3.S) is similar to the single 
curvature bending case, when p < 0.5, i.e. the buckling 

capacity decreases with decreasing a. But for p > 0.5, the 

buckling capacity reduces with increasing value of a. The 
Pe^ component of major axis moment reduces the portion of the 
smaller flange which is under compression. As ci increases, 
Pe^ decreases, thus increasing the portion of the smaller 
flange under compressive stresses. This results in reduction 
in the buckling capacity. 

(ix) Figure 3.9 shows the behaviour of tie-beams under moment 

gradient. For single curvature bending {ft ^ O) , the 

behaviour is similar to the beam case (Fig. 3.7, ft < 0} 

except for very low value of a, where the buckling capacity 
for p > 0.5 reduces with increasing p value. The reason is 
that, for p > 0.5, component increases the effective 

moment and for low values of a it over shadows the Wagner 
effect. 

(x) For double curvature bending of tie-beams {ft = +1.0, Fig. 

3.9), a doubly symmetric section has the highest buckling 
resistance. For all values of p, the buckling capacity 
increases with increasing a. 

(xi) Figure 3.10 shows the influence of K on buckling capacity. 
It is interesting to note that for prismatic member (a = 1.0, 
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Fig. 3.10a), with very high value of p and with large anial 
compression, virtually there is no buckling resistance for K 
> 1.0. However, it is not so for the case of tapered member 
with a = 0.25. As mentioned earlier (Refer Fig. 3.5b, 
sagging moment case), this is because of the presence of Pe_, 

"Z. 

which reduces the portion of smaller flange under 
compression. 

3.3 TfcANSVERSE LOAD CASES 

Results for transverse load cases are presented in Figs. 3.11 

- 3.15. the stability criteria for tapering monosymmetr ic 

members under central point load (assumed at shear centre) is 
given in Fig. 3.11 for K = 1.0, ot = 0.25, 1.0 and p = 0.1, 0.5, 

0.9. Figure 3.12 shows the influence of a. and p on the buckling 
capacity with K = 1.0 and X = -1.0, 0.0, 0.5. The influence of 
load height below the shear centre for K = 1.0, a = 0.25, 1.0 and 
p = 0.5 is shown in Fig. 3.13. the effect of the positioning of 

the single point load along the span can be seen in Fig. 3.14 for 

p = 0.1, 0.5, 0.9; X = -1.0, 0.0, 0.5; K = 1.0 and a = 1.0. 

Similar effect in the case of two point loads is shown in Fig. 

3.15. From these results following points can be noted! 

(i) It is seen from Fig. 3.11 that for X > 0.25, the buckling 

curves are similar to those in Fig. 3.3 for the uniform 
sagging moment case, that is, k increases with decreasing X 
and with increasing p values. However, when the aKial load 
is tensile (X < O), the buckling capacity for p = 0.9 is less 
than that for p = 0.1 and p = 0.5. (This is because of Py 

O "Z 

component which acts in the same direction as the external 
moment when p > 0.5. This component increases the effective 



moment, over shadows the Wagner effect and thus reduces the 
buckling capacity.) A similar trend in the results is 
observed when transverse load is in the form of a uniformly 
distributed load and two point loads. 

(ii> For beams (X = 0.0) the value of k increases with increasing 
value of p, except for p = 0.9 where the buckling capacity 
reduces considerably (Fig. 3.12b). Similar result has also 
been reported earlier for prismatic monosymmet ri c member 
under central point load. The reason is, that for very high 
value of p (T - section), the buckled mode shape for angle of 
twist <p, changes to more complex lower buckled mode shape 
(Wang and Kiti pornchai , 1986b). For beam columns (X = 0.50, 

Fig. 3.1Ec), the behaviour is similar to that for the uniform 
sagging moment case (ft = -1.0, Fig. 3.8). For tie-beams (X = 
—1.0, Fig. 3.12a), the buckling capacity reduces considerably 
for p > 0.5. This is because of predominant effect of ^^ 02 * 
A similar trend in the results is observed when the 
transverse load is in the form of a uniformly distributed 
load and two point loads. 

(iii) It is clear from Fig. 3.13 that the buckling load increases 
as the point of application of the load moves towards the 
bottom flange. The bottom flange loading causes a 

stabilizing torque as the twisting of the cross-section 
occurs during buckling; the top flange loading causes a 
torque which further destabilizes the member. Thus, the 
buckling capacity for bottom flange loading is higher than 
that for shear centre loading and it is least for top flange 
loading. The behaviour is found to be similar for other 


values of p as well. 



(iv) Figure 3.14 shows that for p > 0.5, the buckling capacity of 
a member reduces as the point load moves towards the mid 
span. No significant change in buckling capacity is observed 
for low values of p. The behaviour is found to be similar 
for other values of a as well. 

(v) In case of two point loads and p < 0.5 <Fig. 3.15), the 

buckling capacity of a member increases as the distance 
between the loads decreases (i.e. as the value of a 

increases). For p > 0.5 however, the buckling capacity first 
reduces and then increases. The behaviour is found to be 
similar for other values of ot as well. 


3.4 COMPARISON WITH IS SPECIFICATIONS 

IS:800 (1984) proposes the following formiula for determining 

elastic critical stress f . (in MPa)* 

cb 


f 


cb 


26.5 X 10 ^ 
( 1/ry)^ 
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(3.2) 


where , 


1 

r 

y 



K 


2 


effective length of compression flange; 

radius of gyration of the section about its axis of 

mininium strength; 

mean thickness of the compression flange; 
overall depth of beam; 

respectively the lesser and greater distances from 
the section neutral axis to the extreme fibre; 
a coefficient to allow for the inequality of 
flanges, and depends on w, the ratio of the moment 
of inertia of the compression flange alone to that 
of the sum of moments of inertia of the flanges. 



each calculated about its own axis parallel to the 

Y-Y axis of the girder, at the point of maximum 

bending moment. Value of is given by; 

Kg = 0.5 (au3 - 1) for > 0.5 

= Bui — 1 for Oii < 0.5 

The above code provisions do not take into account the effect of 

web tapering, the loading condition, load height below shear 

centre, the nature of axial force etc. The values of f . 

cb 

calculated as per the code provisions are compared with those 
obtained from present analysis; this comparision is given in 
Tables 3.1 and 3.2 for prismatic beams (ot = 1.0, X = 0.0). The 

details of sections used are listed below: 


No . 

P 

h 

w 

Cmm) 

(mm) 

''t 

(mm) 

®b 

(mm) 

^b 

(mm) 

T 

w 

(mm) 

1. 

0.1 

300.0 

46.0 

38.0 

140.0 

12.5 

7.5 

2. 

0.3 

300.0 

92.0 

19.0 

140.0 

12.5 

7.5 

3. 

0.5 

300.0 

140.0 

12.5 

140.0 

12.5 

7.5 

4. 

0.7 

300.0 

140.0 

12.5 

92.0 

19.0 

7.5 

5. 

0.9 

300.0 

140.0 

12.5 

46.0 

36.0 

7.5 


Note that for all the sections considered fj = 0.5, (B/h^) = 0.467, 

(T/h ) = 0.04 and (T /h ) = 0.025. 
w w w 

Table 3.1 compares f^j^ for prismatic, monosymmet ri c beams 
under uniform sagging moment load case (moment gradient with ft = 
-1.0, p varying from 0.1 to 0.9; 1/ry varying from 100 to 300). 
Table 3.2 shows the comparison for prismatic doubly symmetric 
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beams under various loading cases. A study of above tables 

reveal following important points: 

(i) For p < 0.5, IS code provisions give unsafe value of 

over estimates by 46% to 50% for p = 0.1, and 15% to 18% 

for p = 0.5, depending on the slenderness ratio. 

(ii) For p > 0.5, IS:800 gives conservative values of 
p = 0.9 it under estimates f . by 27% to 40%. 

(iii) For doubly symmietric beams with transverse loads placed on 
the top flange, the value of f^j^ obtained as per code 
specifications is unsafe. 

For the central point load case, Eq. (3.2) over 

estimates f^j^ by 8% to 29%; for the uniformly distributed 
load case the error is 24% to 44%. 

(iv) For doubly symmetric beams with transverse loads placed on 
bottom flange, IS:800 gives conservative values of fj.jj» 

For the central point load case, IS:800 under estimates 
f^j^ by 30% to 45%, while for the uniformly distributed load 
case the error is 12% to 28%. 

Tables 3.3 to 3.10 give values of non dimensional buckling 
capacity k, for doubly symmetric sections under various loading 
cases. The value of a. is either 0.5 or 1.0. The axial load ratio 
varies from -1.0 to +1.0, while K varies from 0.25 to 2.0. Once a 
value of k is picked up from the appropriate table, the f^j^ can be 
obtained as 
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3. 5 SCOPE OF FURTHER RESEARCH 

The present investigation is based on certain simplifying 
assumptions, which are listed in Sec. S.1. Further improvement is 
possible along the following liness 

(i) The buckling capacity of a member may be increased by the 
proper use of a system of lateral bracing. Such systems may 
consist of either discrete braces (e.g. series of cross beams 
connected to the member) or of a continuous lateral restraint 
(e.g. members supporting concrete floor). Most bracing 
arrangements in an I section member can be represented by an 
idealized spring system consisting of an elastic lateral 
brace at a certain distance above the shear centre of the 
cross section and an elastic rotational brace at shear 
centre. Flexural-torsional behaviour can be investigated for 
such bracings placed at various positions along the span. 

(ii) End restraints assumed in the present analysis are! 

- No lateral displacement and moment at supports. 

- Zero twist and zero warping restraint at supports. 

Similar investigation can be carried out for other 
boundary conditions also. 

(iii) The analysis deals with two way tapering members which are 
used essentially for plate girders or gantry girders. 
Sometimes, one way tapering members are used in freae 
construction (e.g. legs of a hinged portal or gable frame). 
The flexural-torsional behaviour of such members can be 
investigated for various types of end restraints. 
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•ig. 3.3(a) Stability criteria for monosymmetric beam-column 
tie-beams under uniform moment, p =-1.0 , d = 1.C 
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Fig.3.3(b) Stability criteria for monosymmetric beam -columns 
tie-beams under uniform moment , p = “1.0, ot =0.25 
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■ig.3.4(b) Stability criteria for monosymmetric beam -columns] 
tie-beams under one end moment, p=0, cL =0.25. 






Axial Load Ratio (A) 


Fig. 3.5 (a) Stability criteria for monosymmetric beam -column 
tie-beams under moment gradient, p = 0.5 , Qt = 1.{ 
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Fig.3.5(b) Stability criteria for monosymmetric beam-column 
tie-beams under moment gradient, p=0.5,c£ =0.2 
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Fig. 3.6(a) Stability criteria for monosymmetric beam-columr 
tie-beams under moment gradient , p = 1.0, c£ = 1.0( 
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Fig. 3.6(b) Stability criteria for monosymmetric beam-columns 
tie-beams under moment gradient, [3 =1.0 , ot =0.25 
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Fig. 3.8 Influence and p on buckling capacity of beam-column 
under moment gradient. 
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Fig. 3.11 Stability cri'eria for monosymmetric beam 
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(a) Tie-beams (b) Beams _ (c) Beam-columns 

( X=-1.0, K=1.0) ( A=0.0, R= 1.0) ( A=0.50, K =1.0) 

Fig, 3, 12 In-fluence of a and P on buckling capacity under central point 
load (SC). 
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Fig. 3.14 Influence of load position along span 
on buckling capacity for single point 
load case (SC) 
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Fig. 3.15 Influence of load position along span 
on buckling capacity for two point 
load case (SC). 







Table 3 


Value Df based on present analysis and 13:800 

for prismatic, doubly symmetric beams 
under various transverse loads. 


Present Analysis 


1/r 13:800 

y 


Central Point Uniformly Distributed 
Load Load 




TF Load 

B L 0 cl d 

TF Load 

BF Load 

•100 

4-19.43 

3?^' . 49 

756.88 

291,45 

585 . 20 

•1S5 

298.81 

239 . 03 

515.67 

212.96 

400.85 

•150 

230.69 

189.86 

382.40 

1 6)B . 36 

298.76 

•175 

187.63 

158.45 

300 . 25 

139 . 92 

235 . 66 

200 

158. 18 

136.66 

245.51 

120.24 

1 93 . 49 

225 

136.62 

120.62 

206 . 85 

105.79 

1 63 . 62 

250 

120.65 

1 08 . 27 

178.29 

94.68 

141 .49 

r--, “7 c 

»v: 

107.98 

98 . 42 

156.44 

85 . 85 

124.51 

300 

97. fj 

90.35 

139.23 

78 . 64 

1 1 1 . 1 0 
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C-.; 

5.20Z-11 

4.38668 

3 . 45385 

sn 

5 . 85E38 

4.80005 

3.63595 

75 

6.69797 

5.36097 

3.91850 

00 

7.68919 

6 . 03595 

4.28189 

25 

8. 784 15 

6.79522 

4.70768 

50 

9 . 95220 

7.61575 

5.18061 

*7 c:._ 

1 1 . 17217 

8.48087 

5.68898 

00 

1 2 . 42966 

9.37886 

6.22413 
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e 3.7 

77 = 0.67 

a — 0 « ; 

jn GJ~ 

fj = 0.5 

p = 0.: 

y 0 0 



X 

0.25 0.50 

0.75 

1 .00 


2.91526 

2 . 29289 

1 . 50778 

-0.27506 

2 . 98446 

2.25128 

1 .35704 

-0,54917 

3.13105 

2.26445 

1 . 23870 

-0.82284 

3.34131 

2.32374 

1 . 14824 

-1.09642 

3.60196 

2 . 42053 

1 .08114 

— 1 . 37 002 

3.90156 

2.54710 

1.03318 

-1.64366 

4.23105 

2.69707 

1 . 00068 

-1.91732 

4.58349 

2 . 86537 

O . 98060 

-2. 19100 
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